Surface properties of fluids of charged platelike colloids 
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Surface properties of mixtures of charged platelike colloids and salt in contact with a charged 
planar wall are studied within density functional theory. The particles are modeled by hard cuboids 
with their edges constrained to be parallel to the Cartesian axes corresponding to the Zwanzig model 
and the charges of the particles are concentrated in their centers. The density functional applied 
is an extension of a recently introduced functional for charged platelike colloids. Analytically and 
numerically calculated bulk and surface phase diagrams exhibit first-order wetting for sufficiently 
small macroion charges and isotropic bulk order as well as first-order drying for sufficiently large 
macroion charges and nematic bulk order. The asymptotic wetting and drying behavior is investi- 
gated by means of effective interface potentials which turn out to be asymptotically the same as for 
a suitable neutral system governed by isotropic nonretarded dispersion forces. Wetting and drying 
points as well as predrying lines and the corresponding critical points have been located numerically. 
A crossover from monotonic to non-monotonic electrostatic potential profiles upon varying the sur- 
face charge density has been observed. Due to the presence of both the Coulomb interactions and 
the hard-core repulsions, the surface potential and the surface charge do not vanish simultaneously, 
i.e., the point of zero charge and the isoelectric point of the surface do not coincide. 
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I. INTRODUCTION 



Clay particles are examples of charged colloidal plates 
with diameters in the range from nanometer to microm- 
eter. Besides their technological importance in applica- 
tions such as agriculture (as soil), construction (in con- 
crete), filling (in cosmetics, rubber, plastic, etc.), coating 
(of paper), and oil-drilling (as rheological fluids for oil re- 
covery), clay suspensions play an important role in basic 
research because they exhibit a complex phase behav- 
ior in terms of packing and orientation of the platelike 

particles [JIlllSHSSSSQilini. Exam P les 
include sol-gel transitions, flocculation, aging, and even 
liquid crystalline phases depending on numerous param- 
eters such as particle size, shape, charge, polydispersity, 
ionic strength, solvent, temperature, and external fields. 
The challenge is to understand the relation between these 
parameters and the phase behavior in order to control the 
properties of such fluids. 

Most studies in this direction have been devoted to 
bulk properties. Some experimental progress has been 
made due to advances in the synthesis of well defined 
and experimentally convenient model systems [l2l Il3| . 
For a recent review of these investigations see, e.g., 
Ref. 14]. Computer simulations for platelike par- 
ticles in bulk systems have been performed in recent 
years 13 111 13 13 ES E3, facing particular 
challenges as compared with systems consisting of rod- 
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like particles. Most theoretical work has been concen- 
trated on the inve stig ation of models of hard platelets 
|22l [231 1241 125L l2a l27l and anisotropically charged discs 
|28l 129 . l30L I31L I32I l33l l34l I35L l3q in spatially homoge- 
neous configurations. 

In contrast to the meanwhile well studied case of bulk 
suspensions of platelike particles, only a few theoreti- 
cal investigations of such systems with spatial inhomo- 
geneities such as interfaces, surfaces, or due to external 
fields have been undertaken [3I I3I liol lil lil li^ 
because of the technical difficulties implied by reduced 
translational symmetry. On the other hand, understand- 
ing the influence of surfaces is of particular importance 
for applications because, e.g., omnipresent walls have a 
strong influence on this sort of material and electrodes are 
a common means to manipulate fluids of charged parti- 
cles. A first important theoretical step for understanding 
such surface properties requires to determine the corre- 
sponding surface phase diagrams. Whereas surface phase 
transitions for simple fluid models have been investigated 
for decades 0, 0, litl l47j , no such studies are available 
for complex fluids involving charges and anisotropically 
shaped particles. 

Here we present the investigation of a density func- 
tional theory for a fluid of charged platelike colloids and 
salt in the presence of a charged hard wall. Due to the 
expected spatial inhomogeneities of the fluid close to the 
surface, density functional theory |4ct is the method of 
choice for addressing these issues |49| . Recently, we have 
proposed a density functional for fluids of charged plate- 
like colloids which also takes long-ranged Coulomb inter- 
actions into account Here a straightforward exten- 
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sion of that functional is presented in order to capture 
the influence of a charged hard wall. 

As introduced in Ref. [i^. the fluid is modeled by 
a mixture of colloidal rectangular plates and cubelike 
salt ions with the platelet orientation restricted to three 
mutually perpendicular directions (Zwanzig model |50|h 
The charges are concentrated in the centers of the par- 
ticles. The bulk phase diagram of this model exhibits 
one isotropic and one nematic phase separated by a first- 
order phase transition [i^. The salt concentration at 
coexistence of the isotropic and the nematic phase is 
larger for the isotropic than for the nematic state [ZsT ]. 
This phenomenon is well-known for membrane equilibria 
[5ll where it is commonly referred to as the Donnan 
effect 53] . By bringing an isotropic and a coexisting ne- 
matic bulk system into spatial contact, an electric double 
layer forms at the free interface. For negatively charged 
platelets, a negatively charged layer forms on the nematic 
side of the free interface whereas a positively charged 
layer is found on the isotropic side ■ The correspond- 
ing electrostatic potential profile changes from one side 
of the interface to the other by the so-called Donnan po- 
tential m. 

The present analysis is devoted to the investigation of 
the model fluid studied in Ref. 0] in contact with a hard 
wall which acquired charges upon releasing counterions 
into the fluid. Section [D] presents the model and extends 
the functional used in Ref. 43] in order to incorporate 
a wall potential. In Sec. II I II the resulting bulk and 
surface phase diagrams are discussed. Further details 
concerning the wetting (drying) behavior of a charged 
wall in contact with an isotropic (nematic) bulk system 
are given in Sec. IIVI l(Vj> . Electrostatic potential profiles 
and the electrostatic surface potential are investigated in 
Sec. IVII Section lYlII summarizes the results. 



II. GENERAL FORMALISM 
A. Definitions 

We consider a fluid of charged platelike colloids and 
salt in three-dimensional space confined by a single pla- 
nar structureless charged hard wall. 

The fluid model used here is the same as in Ref. ^3 • It 
is a ternary mixture of charged hard square cuboids (see 
Fig. ^| with their edges parallel to the Cartesian axes 
(Zwanzig model |H3]) and dissolved in a dielectric sol- 
vent (e.g., water) with dielectric constant e. The solvent 
is treated as a continuum. For simplicity, the charges 
are fixed, monodisperse, and concentrated at the cen- 
ters of the particles. The particles of the first compo- 
nent represent the macroions M of size Dm x Dm x Lm, 
Dm 7^ Lm, and charge Qm- Within the Zwanzig approx- 
imation, macroions can take three different orientations, 
denoted as M x , M y , or M z corresponding to whether the 
-LM-edges are parallel to the x-, y-, or z-axis, respectively 
(see Fig. ^| . The second component consists of salt ions 




FIG. 1: Macroions M are square cuboids of size Dm X Dm x 
Lm, Dm 7^ Lm with charge Qm whereas salt ions S and 
counterions C are cubes of side length Ds with charges Qs 
and —Qs, respectively. The pointlike charges (•) are localized 
at the centers of the cuboids. The macroions can adopt three 
possible orientations M x , M y , and M z corresponding to the 
Z/M-edges being parallel to the x-, y-, and z-axis, respectively. 

S modeled as cubes of side length Ds and charge Qs (see 
Fig.nj. The salt ion charge Qs is chosen to have the same 
sign as the macroion charge Qm if the latter does not van- 
ish. Finally, the third component consists of counterions 
C guaranteeing overall charge neutrality. They are also 
described by cubes with the same side length Ds but op- 
posite charge Qc ■= —Qs (see Fig.^l. Due to the choice 
of the same sign for Qm and Qs, ions of component C 
are counterions in the usual sense for both the macroions 
M and the salt ions S. 

The wall, which confines the fluid from one side, is 
modeled as impenetrable for the fluid particles. More- 
over, it carries a fixed homogeneously distributed surface 
charge of areal density a which, depending on the sign of 
a, can be assumed to have been generated by releasing 
an appropriate amount of ions of either component S or 
component C into the fluid. 

Qi(r), i £ {M x , M yi M Z ,S, C} denotes the number den- 
sity at point r of the centers of macroions with orienta- 
tion M x , VtZ , salt ions, and counterions, respectively. Note 
that the position r 6 V C R 3 with V denoting the system 
volume is a continuous variable in contrast to the orienta- 
tion of macroions, which varies within a discrete set. We 
introduce the abbreviation g(r) := {qm x {y), ■ ■ ■ , £>c( r ))- 

The system under consideration is coupled to two 
particle reservoirs: One supplies neutralized macroions 
(chemical formula C k M, k := > 0) and the other 
neutral salt (chemical formula CS); ^c k M an d Mcs 
denote the corresponding chemical potentials. These 
molecules dissociate upon entering the system: 

C k M — > kC Qc +M Qm , 
CS — ► C Qc +S Qs . (I) 
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These equilibrium chemical reactions lead to the follow- 
ing relations between the reservoir chemical potentials 
(A*C fc M and ^cs) and the particle chemical potentials 
(fii,i € {M x , My, M z , S, C}, ii Mx = VM y = MmJ: 



k^c + Mm, 
Mc + Ms- 



(2) 



B. Density- functional theory 



Configurations of the system are characterized by the 
set of number density profiles g. The equilibrium states 
minimize the grand canonical density functional |48t l49l 


= E/^ ftW(ln(ft(r)) - 1 - mI + K(r)) 



+F cx [ e ] + C/ solf ( f j), 



(3) 



where Vi is the wall (or substrate) potential exerted on 
particles of type i, F c * is the free energy in excess over 
the ideal gas contribution, and U scli (a) is the self energy 
of the charged wall. The latter is independent of g and 
hence merely shifts f2 by a constant value. It has been 
introduced for later convenience. Finally, the reduced 
particle chemical potentials /i* := /i; — In (A|j with the 
thermal de Broglie wavelength Aj (Am x = A^ = Am z ) 
for particles of type i have been introduced. With the 
reduced reservoir chemical potentials 

Vc h M ■= Vc k M - (in (A 3 Mx y J + fcln (A&)) 
Mcs := Mcs - ( In (A|) + In (A 3 C )) , (4) 
Eq. (J2J takes the form 



Mcs — Mc + Ms- 



(5) 



For given reservoir chemical potentials Hc k M an( i McS' 
the particle chemical potentials fi* are fixed by Eq. (JjjJ 
and the constraint of global charge neutrality: 



fd 3 r ^Q; ft (r) + L4|<7 = 



(G) 



where \A\ is the area of the charged surface. 

The Euler-Lagrange equations of the minimization 
problem read 



sn 



= ln( ft (r))- M t + y i (r)-c i (r)=0 (7) 



with the one-particle direct correlation function 

5F CX 

Ci(r) 



(8) 



Equation Q has to be solved under suitable boundary 
conditions in accordance with the bulk states for given 
reservoir chemical potentials Mc fc M an( l Mcs- These bulk 
phases have already been determined in Ref. |43j . 



C. Excess free energy 

As introduced in Ref. the interactions between 

charges are approximated by 



QiQj 



r — r 



(9) 



where the usual Euclidean norm ||r|| 2 = \J x 2 + y 2 + z 2 is 
replaced by the supremum norm HrH^ = max(|x|, \y\, \z\) 
which is the adapted norm in the context of a Zwanzig 
model for cuboids. 

Employing a Debye charging process, F cx can be ex- 
pressed as |43 



pex _ pe 



pcx,c pcx,c 
el corr ' 



(10) 



The first term in Eq. (|10|) is the fundamental measure 
functional introduced by Cuesta and Martinez- Raton J^, 



F c *< h [g] := J d 3 r $(n(r)) 



(11) 



with the weighted densities 

n a (r) = E / d3r ' ^ r - r >'( r ') ( 12 ) 



for a S {0, lx, ly, lz, 2x, 2y, 2z, 3} and the excess free 
energy density 



(n) = -n ln(l-n 3 )+ - 



\qn-2q n 2x n2yn 2z 



q£{x,y,z} 



- n 3 (1 - n 3 ) 2 



(13) 



It describes the free energy due to the hard-core inter- 
actions Uh(r — r') between two particles of type i and 
j at positions r and r', respectively. U^(r — r') is infi- 
nite if the hard cores of these particles overlap and zero 
otherwise. 

The electrostatic part 

Krie] ■■= ~ E / dV / dV QiWQjVMji* O 
y y y 

accounts for the particle-particle Coulomb interaction 
within the random phase approximation (RPA) whereas 
the correlation part 
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F^M := \ Y,j A " r /dV ft (r)^(r')^(r-r')(exp(-^.(r-r'))-l 

+ exp(-t/£(r - r'))G y -(«(r,r'; [g]), ||r - r'HJ) (15) 



with the function 



Gij(K, s) 



drj min(l,Uij(s)r)exp(—y/rjKs)) (16) 



and the screening factor 



using 



«(r,r';[ g l):=i(«(r;[£|)+«(r';|^)) (17) 



k(t; [£]) := ^Ql^sW + gcW) (18) 



represents corrections beyond RPA. 



Since (k, s — > oo) 



-UQ&jK^s- 5 (see Ref. 



|43j1. the contribution %-(s) := U^{s)Gij{n, s) in Eq. 
(|15|l beyond RPA due to the Coulomb interactions has 
the sewne s -6 asymptotics and sign as the contribution 
within RPA of nonretarded dispersion forces. The am- 
plitudes of nonretarded van der Waals forces for ty pica l 
Hamaker constants of clay (e.g., 6.6 fc^T for AI2O3 pt\ \ 
are two orders of magnitude smaller than those of My (s) 
for \Qm\ = e. Moreover, for long distances (s — ^ 00) 
dispersion forces cross over to an s -7 asymptotics due to 
retardation |46| . Dispersion forces are negligible in this 
retardation limit compared to Vij(s). For these reasons 
we have not included dispersion forces into our model. 



D. Wall potential 

Since the thermodynamic limit of globally charge neu- 
tral systems of hard particles with Coulomb interactions 
exists H2,|5|j, i.e., the bulk free energy density exists and 
it depends neither on the shape nor on the nature of the 
boundaries of the system volume V, the system volumes 
V(L) := A{L) x [0~, L] of size 2L x 2L x L are considered 
in the limit L — > 00. Here A(L) is a square in the x-y 
plane of side length 2L, i.e., |^4(L)| = 4L 2 , with (lateral) 
periodic boundary conditions. Furthermore, the charged 
hard surface is chosen as the intersection of V(L) with 
the plane z = 0. The notation CP in the above defini- 
tion of V(L) is used in order to emphasize the fact that 
the charged wall at z = belongs entirely to the system 
volume. 

The wall potential Vi for particles of type i can be 
decomposed as V- 1 + Vf with the hard wall contribution 



V» := 




(19) 



where Sf is the z-extension of particles of type i, and the 
contribution due to the surface charge 



Vf(r) := J d 



3 , Q t o-5{z') 



(20) 



v 

Similarly, the wall self energy is given by 

/-self 



|: 1/dV /dV ■ (21) 



V V 

Introducing the local charge density 

i 

and the electrostatic potential 



one readily finds 



(22) 



(23) 



v 

d 3 r 



3 , <?<W(r') 



r r' 



d 3 r e °(r)^(r). 



(24) 



According to Ref. |43|], the expression for the electro- 
static potential reduces to 



tp(z) = -4 J dz' q q {z')\z- z'\, 

0" 



(25) 



due to the lateral translational symmetry. By differen- 
tiating twice one can show that i/j fulfills the Poisson 
equation 

ip" = -8g Q . (26) 

Combined with Eqs. IjlUI) and 1241) . the density func- 
tional in Eq. © takes the final form 

L 

%] = \ A \ E y*d*ft(z)(M«(*))-l-Mi 



+^) + ^M) + ^M(0) 
+^ X >>]+F C C ^] 



(27) 
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which has to be minimized under the constraint of global 
charge neutrality 

L 

J dz g Q (z) = 0. (28) 
o~ 

The corresponding Euler-Lagrange equations read 

ln(fc(*)) - + V>{x) + QMz) - 4(z) - c c COIIA (z) = 0. 

(29) 



0(z- 5 ). Note that on the right-hand side of Eq. 

for z — ► oo the leading and the next to leading order 

contributions to the density profiles g cancel. 

Interestingly the asymptotic decay proportional to z~ 3 
of Agi(z) within the current model for charged particles 
equals the corresponding one for systems governed by 
isotropic nonretarded van der Waals forces in Ref. |60| . 
This behavior is in sharp contrast to the results ob- 
tained within multicomponent Poisson-Boltzmann the- 
ories which give rise to exponentially decaying density 
profiles |6lj . 



E. Asymptotic behavior 



III. BULK AND SURFACE PHASE DIAGRAM 



Expanding the Euler-Lagrange equations (Eq. f° r 
z — > oo around the bulk values provides the number den- 
sity deviations Ag{z) := g{z) — g(oo) with 

A ft (z^oo) = &(oo)(- AV^{z)-Q % A^{z) (30) 
+ A$(z) + Ac c COIl!i (z)). 

The electrostatic potential difference Aijj depends on Ag 
via the local charge density g® (see Eqs. I)22|l and J2EJ)- 
However, this dependence will not be used for the deriva- 
tion of, c.f., Eq. Inserting Eqs. ijTTJ) and (fT5|> into 
Eq. JSJ leads to 

Acfc-xx^-^A^zJAy (31) 



and 



Ac c COI ^(z -> oo) = J2 (*«(*) ~ &Qi{z)Bn), (32) 



Upon solving the Euler-Lagrange equations (Eq. 1|29|0 
for mixtures of platelike macroions and monovalent salt 
in contact with a charged hard wall, one finds the equilib- 
rium state of this system. Here the equilibrium states are 
calculated as function of the chemical potential Mc W 
the macroion charge Qm, and the surface charge den- 
sity a. The remaining parameters of the model have 
been fixed as (see Fig. [TJ D M = 20£ B « 14 nm, 
L M = i B « 072 nm, D s = £ B « 0.72 nm, \Q S \ = e, 
and salt density gs = 0.1 mM in the bulk. As discussed 
in Ref. [2iJ, the current model yields reasonable results 
only in the range \Qm\ e, i.e., close to the point of zero 
charge of the colloidal platelets. 

It turns out that the equilibrium states exhibit ro- 
tational symmetry with respect to the wall normal (z- 
axis), i.e., gu x = QMy Therefore, the local equilibrium 
structure of the macroions is captured completely by the 
two densities qm x and qm z - Equivalently, the (total) 
macroion density qm := 2qm x + Qm z and the nematic 
order parameter along the z-axis 



where Uj(z) := -24Q|<3|^(oo)k(oo) 4 z 3 . The spa- 
tially constant 5 x 5-matrices Aij and Bij depend only 
on model parameters (particle sizes and charges) and 
bulk quantities (see Appendix A). By defining Mjj := 
%/f?i(oo) + + and since V^(z) — for z > 
imaxj(S'f) (Eq. (O), Eq. is equivalent to 



MijAg^z -> oo) = -QiAip(z) + ^ ty (2 



(33) 



Firstly, Eq. (|33l) implies an asymptotic decay of Agi(z) 
not faster than proportional to z~ 3 ; otherwise Atp(z — > 
) = Kz~ 3 with some spatially constant amplitude K 
would fulfill the relation K+2AQ i k{oo)- i J2j Q"j9j(°°) = 
simultaneously for all i which is impossible because 
the second term depends on i whereas the first does not. 
Moreover, it can be shown that the deviations Agi(z) 
exhibit a purely algebraic decay in leading order. Finally, 
one is led to the conclusion that the asymptotic Euler- 
Lagrange equations (|33|l in conjunction with the Poisson 
equation (12(it leads to the properties Agi(z — > 00) = 



0(z~ 3 ), A^{z 



3G 



) = 0(z- 3 ), and qV(z -> 00) 



sm 



3 QM* 
2 gn 



1 



(34) 



may serve to identify various structures. The definition 
of sm coincides with the well-known scalar liquid-crystal 
order parameter S — (P 2 (cos-d)) = |((cos$) 2 ) — \ for the 
special case of a Zwanzig model, within which the only 
possible macroion orientations are parallel [M z , cost? = 
1) and perpendicular (M x<y , cosfl = 0) to the z-axis. 
Therefore, structures with sm — and sm > are called 
isotropic and nematic, respectively. 

Figure [3 displays the bulk and surface phase diagram 
for the surface charge density a — —QAe/D\ [ in terms 
of the macroion charge Qm and the chemical poten- 
tial difference A/j^ M :— Hc k M ~ McIm- The solid line 
{A[Xq m = 0) denotes the states of bulk coexistence be- 
tween the isotropic (J) phase and the nematic (N) phase, 
corresponding to the chemical potential Mc/Jm a ^ coexis- 
tence. The bulk equilibrium states for A/j,Q kM < and 
A[Iq M > are isotropic and nematic, respectively. 

At isotropic-nematic bulk coexistence with isotropic 
boundary conditions in the bulk (AfiQ M = 0~), two 
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FIG. 2: Bulk and surface phase diagram of mixtures |54B of 
platelike macroions (Dm = 20£b,Lm = £b) and monovalent 
salt (Ds = £b> |Qs| = e) for salt density qs = 0.1 mM in 
contact with a charged hard wall of surface charge density 
a — — 0.4e/.D M in terms of macroion charge Qm and chemi- 
cal potential difference Afi£ kM . The bulk equilibrium states 
for A/j,Q kM < and A/i^ M > are isotropic (7) and ne- 
matic (N), respectively. Isotropic-nematic bulk coexistence 
corresponds to A(Xq m = (solid line). For isotropic bound- 
ary conditions in the bulk (z — > oo) at isotropic-nematic co- 
existence (Afic ,m = 0~), two first-order wetting transition 



(at Q M 



points W~ (at QZ e [-2 x 10~ 3 e,0]) and W H 
[O.le, 0.132e]) have been found. The tolerance intervals are 
indicated by frames of corresponding widths. Complete wet- 
ting by the nematic phase occurs for Qm £ (Qm~ i Qm ) 
upon approaching coexistence from the isotropic side. For 
nematic boundary conditions in the bulk at isotropic-nematic 
coexistence (A/x<^ fcM = + ), a first-order drying transition 

point D~ (at Q M — — 0.35e) has been found. The accom- 
panying predrying line (dashed line) terminates at a critical 

point C£ (at Ql° = -0.87e,A/^ M = 0.049). Complete 
drying by the isotropic phase occurs upon approaching coex- 
istence from the nematic side for Qm < Qm ■ This implies 
that for Q M < Qm < Qm or QYi < Qm < e there is nei- 
ther complete wetting (by the nematic phase) nor complete 
drying (by the isotropic phase). A second drying transition 
point D + appears (not shown) for Qm sufficiently large and 
a sufficiently small so that there is reentrance of complete 
drying for large positive values of Qm- 



first-order wetting transition points W~ and W + have 
been found which could be located only within the inter- 
vals of the coexistence line indicated by frames of corre- 
sponding widths: Whereas W~ (Qff € [-2 x 10~ 3 e,0]) 
is known rather precisely, there remains some uncer- 
tainty with respect to the location of W + (QYi £ 
[O.le, 0.132e]). The corresponding prewetting lines are 



so close to the coexistence line such that they could not 
be resolved numerically. 

Complete wetting by the nematic phase occurs for ap- 
proaching isotropic-nematic bulk coexistence from the 
isotropic side for Qm <= (<2m i )■ The phenomenon 
that only partial wetting is found for sufficiently large 
macroion charges \Qm\ may be qualitatively understood 
as follows: The macroion number density profiles QM m , 
qm v , and qm z close to the surface are influenced by the 
hard-core interactions — which give rise to a preference 
of nematic order close to the wall — , the macroion- 
substrate Coulomb interactions proportional to |Qm|> 
and the macroion-macroion Coulomb repulsion propor- 
tional to \Qm \ 2 ■ The latter dominates for large macroion 
charges \Qm\ leading to a depression of the values of 
the macroion number densities near the wall which in 
turn prevents the growth of a nematic film. For small 
macroion charges |Qm|> Qm x , 0m v , and qm z near the 
surface are determined by the balance between the hard- 
core interactions and the macroion-substrate interactions 
which lead to complete wetting for attractive walls. Fur- 
ther details of the wetting behavior will be given in Sec. 

im 

In Fig. |21 a first-order drying transition point D~ at 
Q M — — 0.35e is found for isotropic-nematic bulk co- 
existence with nematic boundary conditions in the bulk 
(A/j,Q kAI = + ). The first-order character of the dry- 
ing transition at state point D implies the existence of 
a predrying line (dashed line in Fig. [2J , along which the 
excess adsorption of macroions 



M 



dz (qm(z) - Qm(oo)), 



(35) 



which is proportional to the thickness of the emerging 
film, exhibits a finite discontinuity. The predrying line is 
expected to meet the isotropic-nematic bulk coexistence 
line tangentially [62|. Our numerical data are consistent 
with this behavior but not definitely conclusive due to 
the numerical difficulties arising from the large thickness 
of drying films very close to isotropic-nematic bulk coex- 
istence. The predrying line terminates at a critical point 

C D located at state point (Q A f — — 0.87e, Afi CkM = 
0.049). Complete drying by the isotropic phase occurs 
for Qm < Qm u P on approaching isotropic-nematic bulk 
coexistence from the nematic side. As for the wetting sce- 
nario discussed above, the macroion-macroion repulsion 
will prevail over the macroion-surface interactions for a 
sufficiently large macroion charges |Qm|- Therefore, for 
large |Qm|, the formation of an isotropic film is initi- 
ated by the strongly depressed macroion number density 
close to the surface. Hence one expects complete drying 
for sufficiently large macroion charges |Qjw|. For a sur- 
face charge density a = —0Ae/D M (see Fig. |5J, partial 

drying is found for Qm € (Q M > e L i- e -; the expected sec- 
ond drying transition point D + is located in the range 
Qm > e. For smaller surface charge densities a, indeed 
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two drying transition points have been found within the 
range Qm S [— e , e ]- More details of the first-order dry- 
ing transition and the corresponding predrying line can 
be found in Sec. 

Within the intervals Qm € [Q M , QYi ] and Qm € 

■)W+ r>r>+i 



[Qm > Q 

[Qm^ tQa^] on ly partial wetting and drying occurs. In 



order to validate the topology of the bulk and surface 
phase diagram in Fig. [5J we investigated a modified ver- 
sion of the density functional (Eq. in which the cor- 
relation term F c e ^ r c [gi] is omitted. An asymptotic analy- 
sis analogous to Sec. Ill El leads to exponentially decaying 
profiles of Ag, Aip, and g®. The corresponding phase di- 
agram is qualitatively the same as in Fig. EI In particular, 
there are also first-order wetting and drying transition 
points which are separated by intervals of only partial 
wetting and drying. Due to the exponentially decaying 
electrostatic potential, this modified model is similar to a 
recently investigated model of hard rods interacting with 
an exponentially decaying wall potential 63] . Whereas 
the authors of Ref. [63j ] relied entirely on a numerical 
approach, which provided them only with evidences of 
first-order wetting transitions, we shall show analytically 
that the wetting and drying transitions shown in Fig. [21 
are of first order (Sees. IIVI and IV|) . 



IV. WETTING 

In this section, we discuss the wetting behavior of 
the model fluid of platelike macroions and salt in con- 
tact with a charged hard wall. For a wetting scenario, 
the boundary conditions imposed on the solutions of the 
Euler-Lagrange equations (Eq. (j29(l ^ require the isotropic 
bulk structure far from the wall. Isotropic-nematic co- 
existence with such boundary conditions is denoted as 

In principle, the asymptotic wetting behavior can be 
completely inferred from the effective interface poten- 
tial il cS (C) := {^[g ( ] - where n b denotes the 
bulk contribution to the grand potential and the density 
profiles g^ are the solutions of the Euler-Lagrange equa- 
tions (Eq. 1|29|0 under the constraint of a prescribed film 
thickness £ characterized by the position of the isotropic- 
nematic interface |46j |. Since exact solutions are out of 
reach for the present model, we have determined f2 (£) 
approximately by considering the following subspace of 
density profiles |60| : 



Qi,c{z) ■= 



QT n 



(*-0 



,^<C/2 
,z>(/2 



(36) 



where g wn (z) and g nl (z) correspond to the wall- nematic 
and the free nematic-isotropic density profiles with the 
surface and the interface located at z = 0, respectively. 
Therefore, in the interval z 6 (-co, £/2] the trial den- 
sity profile g^ is described by the wall-nematic profile 

whereas in the interval z S (£/2,oo) it is given by the 
free nematic-isotropic profile shifted to position £. Due 



to £>J™(oo) = oo), the discontinuity of g^ at z = (/2 
vanishes in the limit £ — ► oo. The transition regime 
around z = £/2 does not contribute to the leading asymp- 
totic terms of fl cS (( -» oo) |6fT|. 

Substituting the trial density profiles g^ into the den- 
sity functional in Eq. I|27l) and using the asymptotic be- 
havior of g wn and p ni (Sec. Ill Ej) leads to the effective 
interface potential 

n cS (() = lwn + lm - (q m (-oo) - ^(oo))CA Mm 

+a 2 r 2 + a 3 r 3 + o(r 4 ) (37) 

with the amplitudes 02 and 03 being independent of £. 
Since the density profiles decay proportional to z~ 3 to- 
wards the bulk values (Sec. Ill Ef) . the results of Ref. [6(j 
can be used directly: 02 depends only on the particle 
charges and bulk densities whereas 03, in addition, con- 
tains contributions due to the wall-nematic and the free 
nematic-isotropic excess adsorption (see Appendix A). 
The only difference between the density functional in 
Eq. (|27J) and the one investigated in Ref. |^3| is the pres- 
ence of the electrostatic term Eq. (|24f) . It can be shown 
that it merely contributes a term 0((~ 4 ) to fi cff (£) due 
to the more rapid asymptotic decay of the charge density 
qQ(z 00) = 0(z- 5 ). 

By inspection one recognizes the same leading asymp- 
totic decay of f2 cff (£) in Eq. I|37[) as for systems governed 
by isotropic nonretarded dispersion forces [46L |47j • We 
regard this as an surprising result because our model 
(Sec. Ill C|l does not include dispersion forces. The 



asymptotic behavior 



in Eq. l|T7|) is ultimately gen- 



erated by the assumption of a pair distribution function 
<7y which decays asymptotically with a Debye-Hiickel 
(screened Coulomb) form (4^. This assumption con- 
cerning the pair distribution function gij is valid because 
the linearization approximation underlying the Debye- 
Hiickel theory is justified at large distances. 

Minimizing the effective interface potential fl c (() 
in Eq. I|37ll with respect to the interface position 
C leads to the (equilibrium) excess adsorption Tm ~ 

(-A/iCfcM)~^ /" 00 for A Mc fc M / as long as a 2 > 
0. Evaluating the analytic expression for a 2 along the 
isotropic-nematic coexistence line, which depends only 
on bulk quantities (see Ref. [6(| and Appendix A), one 
finds two wetting transition points — corresponding to 
W~ and W + in Fig. El — with complete wetting, i.e., 
ai > 0, in between. 

Figure displays the excess adsorption Tm (Eq. 
I|35jl) as a function of the chemical potential differ- 
ence Afig kM < for fixed macroion charges Qm- For 
An* CkM /* 0, isotropic-nematic bulk coexistence is ap- 
proached along vertical thermodynamic paths in Fig. [21 
In Fig.[3fa) the curves for Qm = and Qm = O.Ie sug- 
gest a divergence of Tm in this limit, i.e., complete wet- 
ting of the surface by a nematic phase occurs for these 
macroion charges. On the other hand, for Qm < — O.Ole 
the excess adsorption remains finite, i.e., partial wetting 
occurs. The complete wetting curves for Qm = and 
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a = -OAe/Dj 



(a) 

Qm = O.le 




0.2 

io~ 4 



IO" 3 10~ 2 
a = -OAe/Dj, 



10- 



0.5 



(b) 



\ Qm = 

V 



\ Q 



0.2 

io- 4 



10~ 3 IO" 2 



10" 



FIG. 3: Excess adsorption Fa/ of platelike macroions along 
vertical thermodynamic paths in Fig. [^characterized by fixed 
macroion charges Q m and parameterized by the chemical po- 
tential difference AfiQ kM < 0, which measures the thermody- 
namic distance from isotropic-nematic bulk coexistence, (a) 
The numerically determined excess adsorption Fm remains 
finite upon 

A^c-fcA/ / for Qm < -O.Ole (see also Fig. gj, 
whereas a divergence is suggested for the macroion charges 
Qm = and Qm — O.le, i.e., there is complete wetting of 
the surface by a nematic film, (b) The comparison of the nu- 
merical solution for Qm = (dashed line) with the analytical 

asymptotic power law behavior Fm ~ (— Anc kM )~ 3 (dash- 
dotted line) (see main text) indicates that the asymptotic 
regime is not yet reached within the numerically accessible 
range of undersaturations. 



Qm = O.le in Fig. |3fa) exhibit no discontinuity in the 
shown range. Thus, the prewetting lines attached to the 



first-order wetting transition points W~ and W + must 
be closer to the isotropic-nematic coexistence line than 
the numerically accessible values of Afi^ M . 

Figure|3Jb) compares the analytically obtained asymp- 
totic behavior for the case Qm = (dash-dotted line) 
with the corresponding numerical solution (dashed line) . 
The differences between these curves indicate that the 
ultimate asymptotic regime is not yet reached within the 
numerically accessible range of undersaturations. 

The order of the wetting transitions at state points 
W~ and W + in Fig. O is determined by f2 cff (£) for 
A/ig, M = 0. Critical wetting occurs for a% = 0, pro- 
vided 03 > 0; if 03 < the wetting transition is of first 
order and does not necessarily occur at the point given 
by a 2 = [Zrl Ir30l | . According to Ref. [6(j and Appendix 
A the analytical expression for 03 contains a contribution 
due to the wall-nematic excess adsorption, which is influ- 
enced by the surface charge density a. Therefore, there 
is the possibility that the order of the wetting transition 
depends on the surface charge density a. For the values 
of er used in the present study, however, we always find 
£13 < 0, i.e., the wetting transitions at W~ and W + are 
of first order. 

Whereas the asymptotical analysis above is reliable 
with respect to the order of the wetting transitions at 
W~ and W + , this is not the case concerning the loca- 
tion of W~ and W + , because the wetting transitions are 
of first order |4fj| |. Therefore, we have to use numerical 
methods. 

In Fig.0]the numerically determined excess adsorption 
at coexistence is shown as a function of the surface charge 
density a. Without loss of generality, only curves for neg- 
ative macroion charges Qm are displayed: since the den- 
sity functional described in Sec. ^is invariant under the 
simultaneous inversion of the signs of all charges (Qm, 
Qs, and ct), the curve for Qm = Qm 1S mapped onto the 
curve for Qm = ~Qm by reflecting it at the axis a = 0. 

As expected, the excess adsorption T M of the 
macroions increases with increasing surface charge den- 
sity a for fixed macroion charge Qm < 0, because the 
surface becomes increasingly attractive (or decreasingly 
repulsive) for the macroions. However, for a fixed sur- 
face charge density a and sufficiently large macroion 
charges |Qm|, Ta/ decreases upon increasing |Qm|, ir- 
respective of the sign of Qm, i-e., even for oQm < 0, 
for which the wall attracts macroions. This depression 
of the macroion number density near the surface occurs 
because the macroion-macroion repulsion dominates the 
macroion-surface interactions (see Sec. Illlll . From Fig.0] 
one can indeed infer that there is partial wetting for suf- 
ficiently large macroion charges |Qm|- 

The variation of the macroion excess adsorption Tm 
upon Qm /* along isotropic-nematic coexistence is 
shown in Fig. [5] Tm is finite for —Qm > 2 x 10~ 3 e. On 
the other hand, Qm = corresponds to hard colloidal 
platelets for which the occurrence of complete wetting 
is well known |39l |. Hence the lower wetting transition 
point W~ in Fig. [3 is located within the range Q% £ 
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FIG. 4: Excess adsorption Tm of platelike macroions at 
isotropic-nematic coexistence with isotropic boundary con- 
ditions in the bulk (A/i^ fcM =0") (see Fig. |5J in terms of 
the surface charge density a. Curves for Qm = Q M and 
Qm — —Q M can be mapped upon each other by reflecting 
them at the axis a — 0. For fixed macroion charge, the ex- 
cess adsorption increases with increasing surface charge den- 
sity. For fixed surface charge density, the excess adsorption 
decreases upon increasing [Qm]> even for attractive walls. 



FIG. 5: Excess adsorption Tm of platelike macroions at 
isotropic-nematic coexistence with isotropic boundary condi- 
tions in the bulk (A^ ^ = CP) for a surface charge density 
a = — 0Ae/D M (see Fig. [21 as function of the macroion charge 
Qm- It is finite within the range —Qm > 2 x I0 _3 e which 
implies the location of the lower wetting transition point W~ 
(Fig. to be within the range Qm e [-2 x I(T 3 e,0]. 



[—2 x 10~ 3 e,0]. In order to locate the upper wetting 
transition point W + one may use the fact that the true 
wetting transitions points can only be located within 
the interval Q M G [-Q* M ,Q* M ], Q*M = 0.1317165(5)e, 
where ±Q*m are the locations of the wetting transition 
points inferred from the asymptotic analysis. This state- 
ment follows from the observation that \Qm\ > Q*m leads 
to a2 < which renders £ — 00 as a local maximum of 
^ eff (C) S3- Together with the numerically found com- 
plete wetting for Qm = O.le one concludes that the loca- 
tion of the upper wetting transition point W + lies within 
the range Q^ + G [O.le, 0.132e]. 

Finally, Fig. H3 displays the increase of nematic film 
thicknesses upon increasing Qm < in terms of the 
macroion density profiles qm- 



V. DRYING 

In this section, we study the fluid composed of platelike 
macroions and salt in contact with a charged hard sur- 
face for nematic boundary conditions at large distances 
from the wall. Isotropic-nematic coexistence with ne- 
matic boundary conditions in the bulk will be denoted 

as &V>c k M = 0+ - 

We have performed an asymptotic analysis of the ef- 
fective interface potential for drying similar to the one 



A»* CkM = 0-,a = -0Ae/D 




FIG. 6: Macroion density profiles qm of a mixture of plate- 
like macroions and salt at isotropic-nematic coexistence with 
isotropic boundary conditions in the bulk (Afic h M = ) f° r 
a surface charge density a = —0Ae/D M (see Fig. |5J. Upon 
decreasing i.e., upon approaching the wetting transition 

point W~ in Fig. [5] a nematiclike film forms at the surface. 



described in Sec. II VI By means of this analysis two first- 
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FIG. 7: Excess adsorption Tm of platelike macroions at 
isotropic-nematic coexistence with nematic boundary con- 
ditions (A/i^ fcM = + ) for surface charge density a = 
—0Ae/D M (see Fig. in terms of the macroion charge Qm- 
The excess adsorption is finite and bounded from below for 

Qm > Qm = — 0.35e, whereas at Qm = Qm it jumps to 
— oo. It has been verified by a comparison of surface tensions 
that the numerical solutions with finite excess adsorption cor- 
respond to equilibrium structures and not only to metastable 

states. The discontinuity of Tm at Q M identifies state point 
D~ in Fig.|!2]as & first-order drying transition point. The lat- 
ter conclusion can also be drawn from an asymptotic analysis 
of the effective interface potential (see main text). 

order drying transition points D~ and D + have been 
found. Due to the first-order character of the drying 
transitions, the loci of D~ and D + have been determined 
numerically (see Sec. IIV|I . Note that the upper drying 
transition point D + is not visible in Fig. [5] because it 
is located at Q M + > e. However, for sufficiently small 
surface charge densities |c| both drying transition points 
have been found within the interval Qm G [— e, e]. In 
contrast to the wetting scenario in Sec. II VI we have been 
able to numerically detect the discontinuity of the excess 
adsorption at the drying transitions. Moreover, the loca- 
tions of the predrying lines in the phase diagrams (e.g., 
Fig. [2| can be numerically determined. 

Figure [7| displays the macroion excess adsorption Tm 
close to a charged hard wall with surface charge density 
a = —0Ae/D M as a function of the macroion charge 

Qm < 0. For — 0.35e = Q M < Qm < 0, the excess ad- 
sorption is finite and bounded from below. It jumps to 
—oo at Qm = Qm ■ The wall-nematic surface tensions of 
the numerical solutions with finite excess adsorption are 
smaller than the sums of the corresponding wall-isotropic 
surface tensions and the isotropic-nematic interfacial ten- 
sions. Therefore, these numerical solutions indeed de- 



FIG. 8: Excess adsorption discontinuity ATm (solid line) of 
a mixture of platelike macroions and salt at the predrying 
line parameterized by the macroion charge Qm for a surface 
charge density a = —0Ae/D 2 M (see Fig. |5J. The excess ad- 

sorption difference vanishes for Qm \ Qm = — 0.87e ac- 
cording to a power law AFa/ ~ (Qm — Q^.f) 13 with the mean 
field critical exponent f3 = | (dashed line). 



scribe equilibrium structures and not only metastable 
states. The discontinuity of Tm corresponds to the oc- 
currence of a first-order drying transition at Q M , which 
is displayed as state point D~ in Fig. [2 

A first-order drying transition is accompanied by a 
predrying line in the surface phase diagram (see the 
dashed line in Fig. |2J), which connects the drying tran- 
sition point D~ with a critical point C£,, located at 

(QaF = -°- 87e , A V*c k M = 0.049). The predrying line 
may be parameterized in terms of, e.g., the macroion 

charge: Afj,^ P kM (Q M ) for Q M G [Q^f,Q M ) denotes the 
chemical potential difference A/x^ fcM for which the ex- 
cess adsorption Tm as a function of Qm and A/z^ fcM 
exhibits a finite discontinuity AFm ■ Figure |S] displays 
this discontinuity as a function of the macroion charge 
Qm (solid line). It vanishes according to a power law 

Ar M - (Qm ~ Qm) with the mean field critical ex- 
ponent /3 = o (dashed line). Beyond mean field theory, 
one expects an exponent j3 = g, corresponding to the 
two-dimensional Ising universality class. 

FigureEldepicts the formation of an isotropic film upon 
approaching isotropic-nematic coexistence for A/j^ M > 
with macroion charge Qm — — 0.6e and surface charge 
density a — —QAe/D M . A finite discontinuity of the film 
thickness upon crossing the predrying line at Aji^ M w 
0.03 can be inferred. 
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FIG. 9: Macroion density profiles qm of a mixture of platelike 
macroioris with Qm = — 0.6e and salt in contact with a surface 
of charge density a — —0Ae/D M upon crossing the predrying 
line at A/j,q M « 0.03 (Fig. [3). At the predrying line, a quasi- 
isotropic film with finite thickness appears at the surface. 



VI. ELECTROSTATIC POTENTIAL 

Whereas the three preceeding sections have been fo- 
cused on the fluid structure, this section addresses the 
electrostatic properties of the surface due to the contact 
with the fluid of charged particles. 

The electrostatic potential profile difference Aip(z) := 
ijj(z) — ip(oo) (see Eq. JUJ) relative to the electro- 
static bulk potential ip(oo) for Qm — — 0.5e at isotropic- 
nematic bulk coexistence with isotropic boundary condi- 
tions in the bulk (Afig M = 0~) is shown in Fig.^| For 
large distances z from the surface, Atp(z) decays propor- 
tional to z~ 3 (see Sec. Ill E() . Upon increasing the sur- 
face charge density a a crossover from monotonic to non- 
monotonic electrostatic potential profiles occurs at a = 0. 
According to Eq. (|2fj|l . the slope of the electrostatic po- 
tential difference at the wall is given by the surface charge 
density: Aip'(0 + ) = —8a. On the other hand, an electric 
double-layer is formed with a negatively charged layer 
on the nematic side and a positively charged layer on 
the isotropic side if a quasi-free interface between the 
isotropic bulk and a nematic film of finite thickness is 
present |43(. Thus Aip approaches its bulk value from 
below, i.e., Aip'(z) > for large distances z from the sur- 
face. Therefore, for a > 0, i.e., A^'(0 + ) < 0, Atp(z) is 
minimal at some finite distance < zq < oo, whereas for 
a < 0, i.e., Atp'(0 + ) > 0, Atp attains its minimal value 
at z = 0. 

In the aforementioned case, a quasi-free isotropic- 
nematic interface is formed at isotropic-nematic coexis- 
tence for a state of partial wetting. Additionally, nematic 



FIG. 10: Electrostatic potential difference profiles Aip relative 
to the electrostatic bulk potential in a mixture of platelike 
macroions with Qm = — 0.5e and salt at isotropic-nematic 
bulk coexistence with isotropic boundary conditions in the 
bulk (Afic k M = 0") (see Fig. EJ. The slope of Aip at the 
surface (z — 0) is governed by the surface charge density a 
whereas the bulk value A^(oo) = is approached propor- 
tional to z~ s from below. Hence upon increasing the surface 
charge, a crossover from monotonic to non-monotonic electro- 
static potential profiles occurs. 



films of finite thickness are present for states slightly be- 
low isotropic-nematic coexistence. Quasi-free isotropic- 
nematic interfaces also occur for nematic boundary con- 
ditions in the bulk close to complete drying. In the latter 
case, the electrostatic potential decreases upon approach- 
ing the bulk because Qm < 0. Moreover, we have found 
that the electrostatic potential is a monotonic function 
of the distance from the surface if no isotropic-nematic 
interfaces form. 

Figure ^] displays the electrostatic surface potential 
Aip(0) as a function of the surface charge density a 
at isotropic-nematic coexistence with isotropic boundary 
conditions in the bulk (A(j,q M — 0~). Curves for Qm = 
Q M and Qm = —Qm can be mapped upon each other 
by reflecting them at the origin (©). As expected, the 
surface potential increases monotonically with the sur- 
face charge. However, the points of zero charge (a = 0) 
do not coincide with the isoelectric points (A^(0) = 0), 
i.e., the curves in Fig. II II do not pass through the origin 
(0), in contrast to Poisson-Boltzmann theories .61]. This 
behavior arises from the hard-core interaction of the par- 
ticles, which leads to a depletion attraction of the larger 
macroions towards the wall on purely entropic grounds. 
Therefore, for a hard and uncharged wall (a — 0), neg- 
atively charged macroions (Qm < 0) are accumulated 
close to the surface leading to a negative electrostatic 
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FIG. 11: Electrostatic surface potential Aip(0) — ^)(0)—tp(oo) 
as a function of the surface charge density a in a mixture of 
platelike macroions and salt at isotropic-nematic coexistence 
with isotropic boundary conditions in the bulk {A^huM = 
0~) (see Fig.[5Jl. All curves are monotonically increasing but 
they do not pass through the origin (0) at (a = 0, Aip(0) = 0) 
if Qm 7^ 0. Thus, the point of zero (surface) charge (a = 0) 
does not coincide with the isoelectric point (A?/>(0) = 0). Note 
that curves for Qm = Qm an d Qm = —Qm can be mapped 
upon each other by reflecting them at the origin. 

surface potential (At/>(0) < 0). 

Our findings show that the sign of the electrostatic sur- 
face potential Atp(0) for vanishing surface charge density 
cr = depends only on the sign of the macroion charge 
and not on the boundary conditions in the bulk or the 
chemical potential difference A/Iq M . 



The effective interface potentials exhibit the same asymp- 
totic behavior as nonretarded van der Waals forces al- 
though our model does not include dispersion forces. The 
origin of this phenomenon can be traced back to the 
application within the density functional of a particle- 
particle pair distribution function which decays asymp- 
totically with a Debye-Hiickel form. 

Complete wetting by the nematic phase occurs in be- 
tween two first-order wetting transition points on the 
isotropic-nematic bulk coexistence curve for isotropic 
boundary conditions in the bulk (Fig. |3J) whereas only 
partial wetting occurs at the remaining state points of 
the isotropic-nematic coexistence line (Fig. 0}. The 
continuous but limited increase of the excess adsorp- 
tion (Fig. |5j) and of the wetting film thickness (Fig. 0) 
upon approaching the wetting transition points is used 
to locate the wetting transition points within the phase 
diagrams. Whereas the first-order character of the 
wetting transitions and the complete wetting behavior 
Tm ~ (— A/x£ M )~3 have been inferred from analyti- 
cal considerations, the location of the wetting transition 
points must be obtained numerically. The corresponding 
prewetting lines could not be resolved. 

First-order drying by the isotropic phase, character- 
ized by a discontinuous divergence of the excess adsorp- 
tion upon crossing the drying transition point, is found at 
isotropic-nematic coexistence for nematic boundary con- 
ditions in the bulk (Fig. 01 . The predrying line, where 
the excess adsorption shows a finite discontinuity (Fig.^JJ, 
terminates at a critical point (Fig. [HJ - 

If quasi-free isotropic-nematic interfaces between the 
bulk and surface films of finite thickness form, a crossover 
is found from monotonic to non-monotonic electrostatic 
potential profiles upon varying the surface charge den- 
sity (Fig. riCj|> . Surface potential and surface charge do 
not vanish simultaneously (Fig. lll|) . i.e., the point of zero 
charge and the isoelectric point of the surface do not co- 
incide due to the presence of both Coulomb interactions 
and hard-core repulsion. 



VII. SUMMARY 

This work has been devoted to the investigation of 
a recently introduced model fluid composed of charged 
platelike colloids and salt (Fig. in contact with a ho- 
mogeneously charged planar wall. The density functional 
introduced in Ref. [43| has been extended by a wall 
potential incorporating a short-ranged hard-core repul- 
sion as well as a long-ranged Coulomb interaction with 
charged fluid particles. By solving the corresponding 
Euler-Lagrange equations for appropriate boundary con- 
ditions, spatially inhomogeneous equilibrium states are 
obtained. 

Bulk and surface phase diagrams (Fig. [21 in terms of 
macroion charge, chemical potential supersaturation, and 
surface charge density have been calculated by means 
of asymptotical analysis of effective interface potentials 
and numerical solutions of the Euler-Lagrange equations. 



APPENDIX A 

The 5 x 5-matrix in Eq. (|31J) is given by 

oo 

A tJ = - I dzc^(z,[g(oo)}), (Al) 

— oo 

where denotes the two-particle direct correlation func- 
tion of the hard core term F cx ' h (Eq. Hill) - ) and g(oo) is 
the set of constant bulk density profiles. 
By defining 

w tJ (z) := J d 2 al^.(a,z)(exp(-[/h(a,z))-l 

A(L) 

+ exp(-l/h(a,z))G y -( K ,||a,z|| 00 )), (A2) 
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which renders Eq. (|15fl in the form 

oo oo 

F^[g] = \\MY.j dz j dz f ei (z)g j (z')v 3ij (z - z'), 



(A3) 



the 5 x 5-matrix in Eq. (|32ll can be expressed as 

oo 

Bij = J dz Wij(z). (A4) 



According to Ref. |60|. the amplitudes a-i and 03 in 
Eq. i|3"7|) are given by 



«2 



45> 



and 



&3 ^^jj(^i,wn <^j.ni) 



(A5) 



(A6) 



with the abbreviation 

T a : = — ^eft-ooXef (-00) - #(«>))■ ( A7 ) 



The quantities 



and 



= dz 1 



gT(*) 



(A8) 



g«(-00) - Qf(00) 

dz(gf(z) - ef(-oo)) + 
dr #(«>)) ] (A9) 



are related to excess adsorptions of the wall-nematic sur- 
face and the free nematic-isotropic interface, respectively. 
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